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Electroweak Sudakov logarithms at high energy, of the form (a/ sin 2 9w) n log m s/M% w , are 
summed using effective theory methods. The corrections are computed to processes involving two 
external particles in the standard model. The results include non-zero particle masses, such as the 
t-quark mass, electroweak mixing effects which lead to unequal W and Z masses, and radiative 
Higgs corrections proportional to the Yukawa couplings. We show that the matching at the scale 
Mw,z has a term at most linear in log s//i 2 to all orders. The effective theory formalism is compared 
with, and extends, previous work based on infrared evolution equations. 
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The Large Hadron Collider (LHC) has a center-of-mass 
energy of ^/s = 14 TeV, and will be able to measure 
collisions with a partonic center-of-mass energy of sev- 
eral TeV, more than an order of magnitude larger than 
the masses of the electroweak gauge bosons. Radiative 
corrections to scattering processes depend on the ratio 
of mass scales, and radiative corrections at high energy 
depend on logarithms of the form log sjM^z- I n high 
energy exclusive processes, radiative corrections are en- 
hanced by two powers of a large logarithm for each order 
in perturbation theory, and the logarithms are often re- 
ferred to as Sudakov (double) logarithms. Electroweak 
Sudakov corrections are not small at LHC energies, since 
a log 2 s/M^ >z /(Aat sin 2 9 W ) ~ 0.15 at y/s = 4 TeV. 
These Sudakov corrections lead to a breakdown of fixed 
order perturbation theory, and need to be summed to all 
orders. 

Electroweak corrections at high energy have double 
logarithms, even for processes which are conventionally 
called inclusive, such as the total e + e~ cross-section at 
large angles, because the colliding particles are not elec- 
troweak gauge singlets [l[ . There are no electroweak sin- 
glet fields in the standard model. A composite particle 
such as the proton, while a color singlet, is not an elec- 
troweak singlet. 

There is an extensive literature on electroweak Su- 
dakov effects @, S 3 1, 0, M ©, M EH, m HI 0, GJ • 

The computations use infrared evolution equations 0], 
based on an analysis of the infrared structure of the per- 
turbation theory amplitude and a factorization theorem 
for the Sudakov form factor [l6j]. These summations 
have been checked against one- loop @, 1, [3 

and two- 
loop [H El El 0, El computations. 

The Sudakov logarithm log(s/M^ z ) can be thought of 
as an infrared logarithm in the electroweak theory, since 
it diverges as Mw,z — * 0. By using an effective field 
theory (EFT), these infrared logarithms in the original 
theory can be converted to ultraviolet logarithms in the 
effective theory, and summed using standard renormal- 
ization group techniques. The effective theory needed is 
soft-collinear effective theory (SCET) [HE1], which has 
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been used to study high energy processes in QCD, and 
to perform Sudakov resummations arising from radiative 
gluon corrections. 

This paper studies high energy electroweak Sudakov 
corrections using SCET, and expands on our previous 
work [19]. In Ref. [l||, we showed how to compute 
log s/My^z corrections to the Sudakov form factor for 
massless fermions using EFT methods. In this paper, 
the results are generalized to massive fermions such as 
the top quark, and include radiative corrections due to 
Higgs exchange. The corrections are computed without 
assuming that the Higgs and electroweak gauge bosons 
are degenerate in mass, as in previous calculations. A 
new feature of EFT matching, the existence of single log- 
arithmic matching corrections [li^ ] , is discussed in detail, 
and proven to be true to all orders in perturbation theory. 
This paper discusses the Sudakov form factor computa- 
tion in detail. The Sudakov form factor is not of direct 
relevance to LHC processes, but it allows us to illustrate 
the EFT method for operators involving two external 
particles. The computations of the Sudakov form factor 
given in this paper can be used to compute electroweak 
corrections to processes relevant for the LHC, such as di- 
jet production, tt production, or squark pair production, 
which involve operators with four external particles. The 
results are given in a future publication [201 ] . and can be 
obtained from the computations given in this paper by 
summing over all pairs of external particles with the ap- 
propriate group theoretic factors. 

The outline of the calculation is given in Sec. [TTJ The 
SCET formalism and the full theory we use for our calcu- 
lations are described in Sec. lIIIl Known results on the ex- 
ponentiation of the Sudakov form factor, and a compari- 
son of the infrared evolution equation formalism with the 
SCET approach is given in Sec. |W] Section W\ discusses 
the calculation of Sudakov corrections for massive gauge 
bosons and massless external particles. Section IVT1 gives 
the proof that there is at most a single logarithm found 
in the matching condition to all orders in perturbation 
theory, and consistency conditions on the matching coef- 
ficients and anomalous dimensions are given in Sec. IVIII 
The extension to massive external particles is given in 
Sec. I Villi for all possible hierarchies of mass-scales, in- 
cluding cases in which particle masses are not widely sep- 
arated, so that multiple scales have to be integrated out 
simultaneously. Massive scalar exchange graphs, relevant 
for Higgs exchange, are computed in Sec. HXl Applica- 
tions of the formalism to electroweak Sudakov corrections 
in the standard model is given in Sec.[X]fbr light quarks, 
the top quark, and leptons. 

Notation: We use a(/i) = a(fj,) / (4ir) , and a,i(n) = 
ai(/i)/(47r) where i = s,2, 1 for the QCD, SU{2) and 
U(l) couplings in the standard model. Hypercharge is 
normalized so that Q = T3 + Y. Logarithms are denoted 
by L A = logA 2 /^ 2 , for A = Q,M,mi,m 2 . C F and T F 
are the Casimir and index for the external particles. We 
use the subscript F for both fermions and scalars, to 
avoid rewriting the same expression twice. 



II. OUTLINE OF CALCULATION 

The physical quantity we study is the Sudakov form 
factor in the Euclidean region, defined as the amplitude 
Fe{Q 2 ) = {P2\0\pi) for the scattering of on-shell par- 
ticles p 2 = mf by an operator 0, with Q 2 = —(p2~ 
Pi) 2 > 0. The timclike Sudakov form factor is given 
by analytic continuation, F(s) — Fe{—s — i0 + ), so that 

log(Q 2 /^ 2 ) -> log(s/M 2 ) " *7T. 

We will compute Fe(Q 2 ) for fermion scattering by 
O = i/jj^tp, ip-ip, V'ct^^Vj scalars scattering by O = 
i^D^cj) - Dftffy, and fer mion to scalar (or vice- 
versa) scattering by O — i/j(f>. All operators are taken 
to be gauge singlets so the incoming and outgoing par- 
ticles have the same gauge quantum numbers, but not 
necessarily the same mass. 

The form factor, Fe(Q 2 ) is computed using a sequence 
of effective theories. For the high energy process con- 
sidered is this paper, there are several widely separated 
scales and we must switch to the relevant theory as we 
move between scales. At scales higher than Q 2 , the the- 
ory is the original gauge theory, referred to as the full 
theory in EFT terminology. The precise theory, and the 
SCET formalism used arc given in Sec. IIIII 

As we move to scales below Q 2 we transition to an 
effective field theory (SCET) where degrees of freedom 
with offshellness on the order of Q 2 are integrated out. 
The full and EFT have the same infrared (IR) physics but 
different ultraviolet (UV) behavior and to ensure that 
the operators in the respective theories have the same 
on-shell matrix elements, we must introduce a matching 
coefficient, exp[C(/x)]. For later convenience, the match- 
ing coefficient is written as an exponential. If the full 
theory is matched onto SCET at \iq then the matching 
coefficient is chosen so that 

(pa|O(M0)bi) - exp[C(/io)] <P2|6(MQ)|Pi) (1) 

where 0{p) is the EFT operator corresponding to the 
full theory operator 0{p). The matching coefficient 
expC(/XQ) is independent of infrared physics, and can 
be computed if perturbation theory is valid at //q. In 
general, a single operator O can match onto a set of op- 
erators Oi in the EFT with the same quantum numbers. 
This occurs, for example, for four-fermion operators in 
the analysis of high-energy parton scattering, and can be 
included by treating all the equations below as matrix 
equations, as is familiar from the well-known analysis of 
operator mixing. The matching coefficient C(/ig) con- 
tains logfiq/Q 2 terms, and there are no large logarithms 
if hq is chosen to be of order Q. We will choose (Xq = Q, 
though any value of order Q is acceptable. Any phys- 
ical observable is independent of the choice for /iq. It 
is conventional to choose c(/x), the coefficient of O in 
the full theory, to equal unity at /x = Q. With this 
choice, which gives the usual normalization for Fe(Q 2 ), 
c(Q) = expC(Q) is the coefficient of O in SCET at 
fi = Q. The evolution of c([x) between scales is given 
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by the renormalization group equation 



fi- 



dfj, 



(2) 



the n-direction, with momentum scaling n ■ p ~ Q, 
n-p ~ QA 2 , p± ~ QA. The ultrasoft field contains gluons 
with all momentum components scaling as QA 2 . 
The EFT fcrmion field satisfies the constraint 



where 'y(fi) is the anomalous dimension of O in the EFT. 

We must repeat this sequence of matching and renor- 
malization group evolution as various energy scales are 
crossed, and more and more degrees of freedom are inte- 
grated out. An advantage of the EFT approach is that it 
divides the full multiscale computation into several sim- 
pler pieces, each of which depends on a single scale. This 
allows one to easily identify which quantities are univer- 
sal, and which ones depend on the specific process. In an 
EFT calculation, the IR divergences in the theory above 
a matching scale must match with the UV divergences in 
the theory below the matching scale. We have checked 
this explicitly for all the computations in this paper. In 
most of the tables, we have given only the finite parts of 
the graphs. 



III. SCET FORMALISM 

SCET is an effective theory that describes energetic 
particles, with energy of order Q, where Q is some large 
scale which characterizes the scattering process. SCET 
contains all the modes of the full theory with invariant 
mass much smaller than Q 2 . The SCET fields and La- 
grangian depend on two null four-vectors n and n, with 
n = (1, n) and n = (1, — n), where n is a unit vector, so 
that n ■ n = 2. In the Sudakov problem, one works in 
the Breit frame, with n chosen to be along the pi direc- 
tion, so that n is along the p\ direction. The momentum 
transfer q has no time component, q° — 0, so that the 
particle is back-scattered. The light-cone components of 
a four- vector p arc defined by p + = n-p, p~ = n-p. In our 
problem, p{ = pxx = v\ = P2± = 0, and Q 2 = p\p2- A 
fermion moving in a direction close to n is described by 
the n-collinear SCET field £„ iP (ir), where p is a label mo- 
mentum, and has components n-p and p± [l7lll8l|. It de- 
scribes particles (on- or off-shell) with energy 2E = n ■ p, 
and p 2 <C Q 2 - The SCET power counting is p~ ~ Q, 



p + ~ QA , p± ~ QA, where A <C 1 is the power count- 
ing parameter used for the EFT expansion. The total 
momentum of the field (, n ,p(x) is p + fc, where k is the 
residual momentum of order QA 2 contained in the Fourier 
transform of x. Note that the label momentum p only 
contributes to the minus and _L components of the total 
momentum. 

The gauge field is represented by several distinct fields 
in the effective theory: n-collincar fields A ntP (x) and n- 
collinear fields An jP (x) with labels, and ultrasoft fields 
A(x) with no label, analogous to the soft and ultrasoft 
fields introduced in NRQCD The n-collinear field 

contains gluons with momentum near the n-direction, 
and momentum scaling n ■ p ~ Q, n-p ~ QA 2 , p± ~ QA, 
and the n-collinear fields contain gluons moving near 



4 Zn,p 



where 



(3) 



(4) 



are projection operators. The leading order fermion La- 
grangian is [H[ 



£",P7T [in- D+ 

2 \ n-p 



(5) 



where iD = id + gA is the ultrasoft covariant derivative, 
and . . . denotes terms involving the collinear gauge field. 
The fermion propagator is 



ifi n ■ p 
2p 2 ' 



(6) 



The effective theory knows about the large momentum 
scale Q through the labels n • P2 and n • p\ on the fields 
t; n ,p 2 an d £,n,pi for the outgoing and incoming particles. 
As a result, SCET anomalous dimensions can depend on 
Q. However, there are no modes in SCET which couple 
n ■ p2 to n ■ pi, so that SCET does not contain modes 
with off-shellness of order Q , which are present in the 
full theory. 

We will also need to introduce SCET fields to describe 
energetic scalar particles, such as the Higgs boson. We 
will use $„ iP as the n-collinear field for a scalar particle 
moving in a direction close to n, analogous to for 
fermions. The field $ n ,p is normalized the same way as 
the full-theory field 0, and produces scalar particles with 
amplitude unity. The scalar kinetic energy term becomes 



D^D»<t> - < p [(n-p) (in ■ D) +p\] $„,. 



(7) 



in the effective theory. It is also convenient to use a 
redefined scalar field, 



<Pn,p = V( n -P) ®n,p, 



in terms of which the kinetic term becomes 



L = <Pli,p \ in - D 



A. 

n • p 



(8) 



(9) 



and has the same normalization as the fermion La- 
grangian Eq. ([5|). The rescaled scalar propagator is now 



1 

p 2 



n ■ p 
p2 



(10) 



0n,p produces scalar particles moving in the n-direction 
with amplitude \J(n ■ p). 



4 




FIG. 1: Graphs in the standard model involving both Yukawa 
and gauge couplings which have no analog in the toy example. 

The theory we consider is a SU{2) spontaneously bro- 
ken gauge theory, with a Higgs in the fundamental repre- 
sentation, where all gauge bosons have a common mass, 
M. This is the theory used in many previous compu- 
tations [J, H, H, 0, [IH, and allows us to compare with 
previous results. It is convenient, as in Ref. [la], to write 
the group theory factors using Cp, Ca, Tp and np, where 
2np is defined to be the number of weak doublets Weyl 
fermions. We will consider this theory with fermionic 
and scalar matter fields in arbitrary gauge representa- 
tions, with the fermions assumed to be vector-like. These 
fields are the external particles in the operators O. We 
will also need to consider graphs which are analogous to 
Higgs exchange graphs in the standard model. For this 
purpose, we will add a gauge singlet scalar field x, which 
couples to the fermions and scalars via gauge-invariant 
interactions, 

iint = —h^,ix4>iipi — hfaiX^lfti, (11) 

h^, i is dimensionless, and h^i has dimensions of mass. 
We will assume that h$,i is independent of Q for power 
counting purposes. In our toy example, x is a gauge sin- 
glet field, and does not break the gauge symmetry. The 
fcrmion masses are independent of the Yukawa couplings 
of x- The toy example Higgs field is a doublet, and breaks 
the gauge symmetry, but does not couple to the matter 
fields. In the standard model, the Higgs field breaks the 
gauge symmetry, and also has Yukawa couplings which 
generate fermion masses. 

The computations are extended to the SU(3) x 
SU(2)l x U(1)y standard model in Sec. [Xj including 
Higgs exchange corrections and unequal gauge boson 
masses. Our results are given to leading order in the EFT 
power counting, i.e. we neglect power corrections of the 
form mf/Q 2 , and M 2 /Q 2 , while retaining all logarithmic 
corrections log m 2 /Q 2 and logM 2 /Q 2 . The gauge boson 
exchange graphs can be obtained from those of the toy 
model. In the standard model, the Higgs field breaks the 
gauge symmetry, and also has Yukawa couplings to the 
chiral fermions. At one-loop, we can obtain the Higgs 
exchange corrections from the ^-exchange graphs in our 



1 This convention for np is used in Ref. [l5l | . Note that the results 
only hold for Ca = 2, since for an SU(N) group with TV > 2, 
a fundamental Higgs does not break the gauge symmetry com- 
pletely. 



toy example. Graphs with Higgs bosons coupling to both 
the fermions and the gauge bosons start at two-loops (see 
Fig. CO). 

IV. EXPONENTIATION 

We start by summarizing some known properties of the 
Sudakov form- factor [22| for the vector current. We will 
see later how the same expressions can be rederived using 
renormalization group methods in SCET. The Euclidean 
form- factor Fp{Q 2 ) has the expansion (L = log(Q 2 /M 2 )) 

F E = 1 

+a(k 12 L 2 + fcnL + fcio) 

+a 2 (fc 24 L 4 + k 23 L 3 + k 22 L 2 + k 21 L + k 20 ) 

+a 3 (fc 36 L 6 + ...)+... , (12) 

with the a n term having powers of L up to L 2 ". In 
the literature, the highest power of L is called the LLp 
term, the next power is called the NLLp term, etc. We 
have included the subscript F (for the form-factor) to 
distinguish it from the renormalization group counting 
described below. 

The series for \ogFp(Q 2 ) takes a simpler form 

logFe = a(fci2L 2 + fciiL + fci ) 

+a 2 (fc 23 L 3 + fc 22 L 2 + fc 21 L + fc 2 o) 
+a 3 (fc 34 L 4 + ...)+..., (13) 

with the a n term having powers of L up to L" +1 , and 
the expansion begins at order a. Note that Eq. (TT3"]) 
implies non-trivial relations among the coefficients k nm 
in Eq. (|12|) . At order n, there are 2n+ 1 coefficients k nm , 
< to < 2n in Eq. (| 1 2[) . but only n + 2 coefficients k nm , 
< to < n + 1 in Eq. (fT3]) . 

The right-hand-side (rhs) of Eq. (fT3|) can be written in 
terms of the LL series L/o(aL) = k\ 2 cA? + k 2 3a 2 L 3 + . . ., 
the NLL series /i(aL) = kna\-+k 22 a 2 L 2 + . . ., the NNLL 
series af 2 (aL) = k w a + k 21 a 2 V. + ... etc. as 

lo gj F B = L/ (aL) + /i(aL) + a/ 2 (aL) + . . . . (14) 

fo and /i begin at order a, and the remaining f n begin 
at order one. 

In this paper, LL, NLL, etc. (with no subscripts) will 
refer to the counting for log Fp . This is also the counting 
appropriate for a renormalization group improved com- 
putation, and is different from the conventional counting 
discussed above. If one looks at the order a 2 terms, for 
example, the conventional counting is that the L 4 term 
is LLp, the L 3 term is NLLp, the L 2 term is N 2 LLp, the 
L term is N 3 LLp, and the L° term is N 4 LLp. Using our 
counting, the terms are given by exponentiating logi^ 
to LL, NLL, N 2 LL, N 2 LL, and N 3 LL, respectively. At 
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higher orders, the mismatch in powers of N between the 
two counting methods increases. 

For precision electroweak studies, the first few orders 
are sufficient. Typical loop corrections are suppressed by 
a /(An). There can be large coefficients in the perturba- 
tion expansion. For example, there are large coefficients 
in the cusp anomalous dimension (see Eqs. (|50I54[) ). In 
this paper we have computed corrections to the Sudakov 
form factor; for dijet production and processes involving 
four-particle operators, the anomalous dimensions are at 
least twice as large as for the Sudakov problem. For 
these reasons, we use the estimate a instead of a/ (in) 
for the size of loop corrections. For QCD, a ~ 0.1, and 
for electroweak corrections, a — > a em / sin 2 6w ~ 0.03. 
logs /Ml ~ 8 for s ~ 4 TeV, so «L ~ 1 for QCD and 
~ 0.2 for electroweak corrections. The NLL series is of 
order ten percent for QCD corrections, and a few per- 
cent for electroweak corrections. The NNLL series is of 
order a percent for QCD, and sub-percent for electroweak 
corrections. 



A. Infrared Evolution Equations 

An expression for the Sudakov form-factor in the limit 
M/Q <C 1 with onshell massless fermions . p 2 = p|_= 0, 
obtained using the evolution equations is jl5l. I22I l23l |24| 



logF £ (Q 2 



ay 

M* A* 2 



: log F Q (a(M)) 
C(a(/*)) +£(a(M)) 



M du' 2 



(15) 



in terms of functions Fq, £, £ and T of the coupling con- 
stant, which have the expansions 

F (a) = l + F (1) a + F Q {2) a 2 + ... , 

r(o) = r« a + r( 2 ) a 2 + ... , 
C(o) = c (1) « + C (2) a 2 + ..., 

£(a) = £«a + £ (2 V + . . . . (16) 

The superscript gives the loop-order of the Feynman 
graphs which contribute. T is known as the cusp anoma- 
lous dimension. 

The gauge coupling constant g satisfies the renormal- 
ization group evolution equation 



T6tt 2 



(16tt 



(16vr 2 ) 



2\3 



(17) 



or equivalently, 



da a I \ 
dfj, = Ma) 



The one-loop coefficient is 

11 4 1 

bo = -g-CU - -T F n F - ~^ T F n s (19) 

where n s is the number of complex scalars, and 2n F is 
the number of fermion weak doublets, in the convention 
of Ref. [H. 

After combining Eqs. p^HT5|) and expanding to order 
a(M) 3 , the form factor takes the following form: 



log F E (Q 2 ) = 

F « + (cW + ^)L + irWi 



a(M) 



~(f) 2 + F„ (2) + (( (2) +^)L 
+ \ (r< 2 ) - 6oC (1 )) L 2 - I (& r«) L 3 | a(Mf 

+ [\ (F (1) ) 3 -^o (2) + ^o (3) + (c (3) + e (3) )L 

+ 1 ( r (3)-6 1 C (1 )-26oC^)L 2 + i(-6 1 r( 1 ) 

a(Mf 



(20) 



2fo (6 C (1) -r( 2 )))L 3 + ^6 2 rWL 4 



A comparison of this expansion with Eq. (|14[) 
shows that /o is determined by r*- 1 ) and 
bo, fi by r( 1 . 2 ),C (1) ,C (1) ,&o,i I and f 2 by 
F^\r^' 2 ' 3 \C (1 ' 2) ,^ 2) ,&o,i,3- In general /„ is deter- 
mined by and terms up to F (n ~ 1] ,T^ n+1 \ (^ n \b n . 
The expression Eq. (|15p is not unique. The identity 



\f ^^^-("M = G(a( Z ))-G(a(y)) (21) 

can be used to show that Eq. (|15p is invariant under the 
transformation 



r(a(/0) - TW)) + dG{a y )] PaW)), 



C(a(/i)) - CKM)) - 2G(a( M )) , 
£(a(M)) -> £(a(M))+2G(a(M)) 



(22) 



As a result, T, £ and £ are not uniquely determined from 
F E (Q 2 ) by Eq. (Tn 



-2a (b a + b x a 2 + b 2 a 3 + . . .) . (18) 



B. Renormalization Group Evolution Equations 

The corresponding expression for F E (Q 2 ) in the EFT 
formalism, as will be derived in Sec. El is 

log F E (Q 2 ) = C(a(Q)) + D (a(M)) 
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+D 1 (o(Af))lo. 



M M 



M 2 



A(a(/i))log^ 



S(a(/x)) 



(23) 



where expC(a(Q)) is the multiplicative matching co- 
efficient at Q 2 , 7(/i) = A(a(^))log(^ 2 /Q 2 ) + B(a(fi)) 
is the SCET anomalous dimension between Q an( i 
M, and expD(a(M)), D(a(M)) = D (a(M)) + 
Di (a(M)) log Q 2 /M 2 is the multiplicative matching co- 
efficient at M. The matching coefficient C and the SCET 
anomalous dimension 7 are independent of physics at the 
low scale M, and so do not depend on the gauge boson 
and Higgs masses. The new feature of the massive gauge 
boson calculation is the existence of a single-log term, 
Di(a(M)) in the matching at M. That there are no 
higher powers of logQ 2 /M 2 in the matching is proved 
to hold to all orders in Sec. IVIl A, B, C, Do,i have loop 
expansions analogous to Eq. ([16)) . The N"LL series for 
log Fe requires A^ n+1 \ B^ n \ £>W, and C^" 1 ). £>W 
contributes only to the a n L term in /„. 
The identity Eq. and 



[Q d M 2 d^' 2 [* d^ 2 Q 2 

Jm 2 i 1 Jm 2 fJ> Jm 2 A 1 A* 

(24) 



can be used to show that Eq. 
redefinitions 



23J) is unchanged by the 



A(o( M )) 
S(a(/x)) 

C(a(Q)) 
D (a(M)) 
£>i(a(M)) 



A(a(A*)) - 
S(a(/x)) 



ggW) 
dH{a^)) 



da 

+2G(o(/x)), 
C(a(Q))+#(a(Q)), 
D (a(M)) - H(a(M)) , 
D 1 (a(M)) + G(a(M)) . 



Pa{a{n)) , 
/?a(a(A<)) 



(25) 



Transformations such as these can arise from a change of 
scheme in the computation of the SCET matching coef- 
ficients and anomalous dimensions. 

We can now demonstrate the equivalence of the Su- 
dakov form-factor in Eq. (123|) and the form-factor given 
in Eq. {TSJl. By taking G{a) = -f(a)/2 in Eq. (22J), 
#(a) = -C(a) and G(a) = -Di(o) in Eq. ||25J|, brings 
Eq. (jT5j) with Eq. (|23|) to a common form, and gives the 
identifications: 

1 mi \\ 1 dD i{ a {^)) i 

- o « M a W) 

2 2 oa 

= r ( a (w) - g fa — AM a O)) > 



C(a(M))+e(«(M)), 



■A(a(/i) 



C(a(M)) + £> (a(M)) = logF (a(M)). 



(26) 



The lhs of Eq. (|26[) is invariant under Eq. (f25j) . and the 
rhs under Eq. (|2"2"]l. The computations of the SCET 
anomalous dimension and the cusp anomalous dimension 
in the literature use the same scheme, so that V = A/2, 
and 



T(a). 
Di(a) = £(a), 



-5(a) 



1 dC(a) 



/3a (a) 



2 9a 
C(a)+A,(a) 



C(o), 
logFo(a) 



(27) 



The expansion of log Fp(Q 2 ) to order a(M) 3 using the 
SCET form is: 



log F E (Q 2 ) = 



a(M) 



+ 



+ 



C (2 



(2) 



(2; 



+ \ (>> + b B^) L 2 - 1 (& A«) L 3 | a(M) 2 

+ J Df ) )L+i( J 4( 3 ) 

26 (s( 2 )+26 C( 1 )))L 2 -l(6aW 



-26^S (1) + 26 A (2 



-6 2 A(!)L 4 
24 



a(M) 3 . (28) 



V. MASSLESS EXTERNAL PARTICLES 

In this section we calculate the form- factor log Fe{Q 2 ) 
for the case Q 2 ^> M 2 3> m 2 , m\. At scales /x > Q we use 
the full theory, and the renormalization group evolution 
of c(/i) is given by 



, dc(/i) 
d/i 



7F(a(A*)) c(fx) 



(29) 



where 7_f(o) = jp^a + Jp a 2 + . . . is the full theory 

anomalous dimension for O. The one-loop values 7^' 
are given in Table HI The general form for Fp given in 
Sec. IIVI is for the vector current, where jp — 0, and 
c(fi > Q) is chosen to be unity. It also holds for the 
other operators with c(fi = Q) = 1 in the full theory. 

The full theory is matched onto SCET at a scale /j, 
of order Q. The effective theory has modes with off- 
shcllness of order Q integrated out, so the matching co- 
efficient depends on log Q 2 / fi 2 , and these logarithms are 
not large if /1 ~ Q. 



7 



o 




C«(/x)/CV 




dM(h)/c f 


11)11} 




~ L Q + "T 


- 2 




4 Id — 6 


-L 2 m - 


h2L M L Q -3L M + 


9 
2 


5,r 2 
6 


V I V 


o 


-L Q + 3L Q + 


,r 2 
6 


8 




— Lm ^ 


h2L M L Q -3L M + 


9 
2 


S^ 2 
6 




2 


-L Q + 4L Q + 


7T 2 

6 


8 


4L Q -6 


-L 2 M - 


h2L M L Q -3L M + 


9 
2 


5^ 2 
6 




-6 


~Lq + Lq + 


,r 2 
6 


2 


4L Q -8 


-L 2 M ^ 


h2L M L Q -4L M + 


7 
2 


5,r 2 
6 







-L Q + 4L Q + 


n 2 
6 


8 


4L Q -8 


-L 2 M - 


h 2I_mLq -4L M + 


7 
2 


5,r 2 
6 


-00, (p'lp 


-3 


-l4 + 2L Q + 


7T 2 
(i 


4 


4L Q -7 


— Lm 


- 2LmLq — |Lm + 


4 - 


6 



TABLE I: One-loop corrections to the Sudakov form-factor, "/f is the full theory anomalous dimension, C is the matching 
coefficient at fj, ~ Q, 71 is the SCET anomalous dimension, and D is the matching coefficient at fi ~ M. 7^ 1 \ C' 1 ', 7J 1 ' and 
D* 1 ' are the coefficients of a = a/(4n) in the one-loop corrections, and Lq = logQ 2 //i 2 , Lm = logM 2 //i 2 . 



The operator O in the full theory matches to the op- 
erator 6 in SCET: 

-» expC( M )[6 1>P2 w„]r[wi^, pl ], 

0V - e X pC( M )[$t iP2 W„][Wi$ R)Pl ], 
z0ty> _> expC( M ) [*J, lPa Wn](*1>i + ^^?2) Al [W^$ fl , ^)1 ] , 
^0 - expC(n){Z niP2 W n }[W^ Pl }, (30) 

where i£>i = V +g(n ■ Ah,,)|, z£> 2 = + ff(n • K,- q )% 
and "P are the SCET label operators introduced in Bauer 
et al. [18| . Collinear gauge invariance requires that, in 
the matching of gauge invariant operators at leading or- 
der in the power counting, the fields occur in the combi- 
nation W^ n<p , W^ n ,p, where W n is a Wilson line con- 
taining n-collinear gauge fields obtained by integrating 
over a path in the n-direction [l8| . 

C(/i) depends on the operator being matched (i.e. the 
C's in Eq. ([50]) have different values, and C can depend 
on r) and, for convenience, we have written the mul- 
tiplicative matching coefficient as exp C ((a) rather than 
C(/i). As is well-known, the matching coefficient can be 
computed as the finite part of the full theory matrix ele- 
ment, evaluated on-shell, with all infrared scales such as 
the gauge boson mass set to zero (see e.g. [H, HI, H3)- 
The full theory graphs to be evaluated at one-loop are 
those in Fig. [21 and when combined with the wavefunc- 
tion and tree- level graphs, give the value of the full the- 
ory matrix element (p2\0\pi). The graphs for the EFT 
vertex correction are shown in Fig. [3J and when com- 
bined with the tree-level and wavefunction graphs, give 
the EFT matrix element (p2\0\pi) . The gauge boson and 
fermion masses are infrared scales and can be set to zero 
in the matching computation thus leading to scaleless in- 
tegrals for the one-loop EFT and wavefunction graphs. 
Since these scaleless integrals are set to zero in dimen- 
sional regularization, the EFT matrix element is equal 
to its tree-level value. The full theory and EFT opera- 
tors O and O are normalized to have the same tree-level 
value, so exp[C(/i)] = (p 2 |C|pi) / (P2|C|pi) troc , i.e. the 



Pi P2 




(a) (b) 



FIG. 2: Graphs contributing to the matching condition 
C(a(Q)). The solid line can be either a fermion or scalar. The 
second graph only exists for the scalar case O = {((p'D^cf) — 



matching condition exp C is given by the on-shell full 
theory matrix element normalized to its tree-level value 
(see e.g. Ref. (2(| for more details). 

The computation of the SCET one loop graphs for O = 
ipj^ip is identical to that for DIS [25|. Particle masses, 
such as the gauge boson mass, are all much smaller than 
Q, and only contribute M 2 /Q 2 power corrections at the 
scale Q, which are being neglected. The one-loop val- 
ues of C(n) for the other cases are computed similarly, 
and are given in Table H where C(p) = a(/i) / '(4tt) 
defines the one-loop correction C' 1 ). There are no large 
logarithms in this matching correction if the matching 
scale /i is chosen to be of order Q. We will choose the 
matching at the high scale to be at fx = Q, and C(/i = Q) 
is given by the third column in Table U with Lq — > 0. 

The renormalization group evolution of c(/^) in the ef- 
fective theory is given by the anomalous dimension of O 
in SCET. The anomalous dimension 71 is used to evolve 
c(/i) from /J, — Q down to the low scale /1 = M . The 
one-loop anomalous dimension is given by the ultraviolet 
countcrterms for the SCET graphs in Fig. [3] (after zero- 
bin subtraction, see Ref [28]). As noted earlier, anoma- 
lous dimensions in SCET can depend on Q. Ultraviolet 
divergences do not depend on the infrared properties of 
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P2 



(a) 



(b) 



P2 



Pi .< 



(c) 



FIG. 3: SCET graphs for the matrix element of O. The dotted 
lines are SCET propagators, and represent either fermions or 
scalars. The upper graphs are the n-collinear and n-collinear 
graphs, and the lower graph is the ultrasoft graph. There are 
also wavefunction graphs. For i(jr D^cf), graph (a) also has a 
contribution where the gauge boson field at <g> arises from the 
covariant derivative. 



the theory, such as a gauge boson mass, so the anoma- 
lous dimension for O = i/jj^iJ) is identical to the DIS 
result . The same argument as that given in Ref. [25[ 
for deep inelastic scatterin g sh ows that 71 is linear in 
log/z 2 /<5 2 to all orders [2al29l]. so 7 is written as 



7i(/f) 



+ B(a(n)), (31) 



which defines A and B. The anomalous dimension has 



the expansion 71 = 7^ 



(!), 



A^a 2 



(2) 2 

7i a 



A^ 



, B = B^a + B^a 2 + ... The computations 
for the other cases are similar, and the results are given 
in Table |U Note that the anomalous dimension depends 
only on the external fields for the operators, and is equal 
for the three fermion operators, and for the two scalar 
operators. The reason is that the EFT anomalous di- 
mension depends on the IR divergence of the full theory 
graph, and the IR divergence is independent of the vertex 
factors. The anomalous dimension for ipcj> is the average 
of the anomalous dimensions for the fermionic and scalar 
operators. 

The next step in the EFT computation is the matching 
condition at the low scale /x ~ M. At this scale, the 
massive gauge boson is integrated out, and one matches 
to an effective theory which is SCET without the massive 
gauge boson. In our toy example, this effective theory 
contains no gauge particles, and is a free theory. There is 
no need to introduce any propagating gauge modes below 
M [3(J. The matching at fj, ~ M is given by evaluating 
the graphs in Fig. [31 and the wavefunction graphs. The 



gauge boson mass can no longer be set to zero, since it is 
of the same order as the matching scale, and the one loop 
SCET graphs are non-zero. The matching computation 
is discussed in detail here for the fermion vector current. 
The other cases are treated similarly. 

One matches the operator c(/j)[£ nj)92 WTi^lW^ft.pJ 
in SCET with gauge particles (the theory above M) 
onto the operator [expD(/i)] c(/i)^ njP2 7 AI ^ ftiPl in SCET 
without gauge particles (the theory below M). 
collinear graph in Fig. [3] gives 



The 7i- 



In 



-ig 2 fi 2e C f 



c{n) 



d d k 



(27r) d k 2 - M 2 



$_ a ji n ■ {p2 - k ) 
2 n 2 (ft - k) 2 ' 

-2ig 2 C Fl ^ 2 



-T- 

d d k 



(27r) d 
n ■ (pa 



k) 



[(P2 - k) 2 + i0+] [-n ■ k + i0+] [k 2 - M 2 + i0+] ' 

(32) 

This integral is divergent, even in 4 — 2e dimensions with 
an off-shellness, unlike the previously studied examples 
where the gauge boson was massless. A related diver- 
gence was encountered by Beneke and Feldman in their 
study of the B — > iriv form-factor. Beneke and Feld- 
man used an analytic regulator [3ll. |32| to evaluate their 
integrals, and we use an extension of their method. A 
similar procedure was used by Jantzen et al. Q in their 
study of two-loop electroweak Sudakov corrections. The 
Pi propagator denominator (jpi — k) 2 in the full theory is 
analytically continued to 



1 



2\<5, 



(Pi - kf 



l+Si 



(33) 



where Vi and 8i are new parameters. The (p2—k) 2 denom- 
inator in Eq. (|32[) arises from the collinear P2 propagator, 
and so gets modified as in Eq. (|33[) . The — n ■ k propaga- 
tor in Eq. (f3"2"| arises from the (pi — k) 2 propagator when 
k becomes n-collinear. In this limit 



1 



2\<5i 



(Pi - kf 



[(n-pi)(-n ■ k)] 



l+Si 



(34) 



We will therefore analytically continue the — n ■ k propa- 
gator in Eq. (|3"2")1 , which arises from the W n Wilson line 
in O using 



1 



5i 



• k 



(-n-k) 



l+Si 



(35) 



where 



v\jp\. We will see below that it is im- 



portant that i>i is related to v{ in this way. Note that 
under boosts, transforms like the minus component 
of a four-vector. With this choice, Eq. (|32f gives 



4 = 



j4_\ 

M 2 ) 



,P2 



y T(e + S 2 ) r(2 - e - <S 2 )r(,S 2 - fr) _ 



T(i + <y 2 ) 



r(2 



Si) 



The regulated value of /„ is given by setting Si = Ti8 and 
taking the limit S — ► first, followed by e — > [3lll32T|. 
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l0g AP 
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(37) 



which is a boost invariant expression, since Vi /p 2 is 
boost invariant. Equation (|37p is valid away from the 
symmetric point n = r 2 . 

The n-collinear graph is given by Eq. (|37p with the 
replacements 5i «-> <5 2 , ^ 2 — ► i/x, — > >^ 2 , P 2 Pi, 



with 



,+ = 



T(e + S 1 )T(2-e-S 1 )T(S 1 -S 2 ) 



r(2- 



(38) 



The parameters v 2 and i/y play the same role as fx in 
the rapidity regularization method of Ref. [28| . 

The ultrasoft graph in Fig. [3] is regulated by the same 
method. The p 2 propagator (p 2 — k) 2 is multipole ex- 
panded in the effective theory, and becomes —p 2 k + , 
where p 2 is a label momentum. Using Eq. (|33[) for the 
fermion propagators, we see that after multipole expan- 
sion, they are regulated in the same way as the Wilson 
line propagators. The ultrasoft graph gives 
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FIG. 4: One-loop correction to the fermion propagator. 
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TABLE II: One-loop gauge boson contribution to on-shell 
wavefunction renormalization. The gauge boson mass is M 
and the particle (fermion or scalar) mass is m. Hf,s are given 
in Appendix [Bl 



n a ■ 



1 1 In M2 

e uv 2 fi 2 



and so contributes a wavefunction correction 



Sz = C F — 
Air 



1 



--L M 



(40) 



(41) 



d d k 



(2Tr) d k 2 - M 2 

u {-"if 1 



[n- { P2 -k)] 1+s * ' [n ■ (pi - fc)] 1+5 i ' 



(39) 



and vanishes on-shell, since p 2 = P\ = 0. 

The total SCET contribution is given by the sum of 
the n-collinear, n-collinear and ultrasoft graphs, as well 
as the wavefunction renormalization correction. The 
collinear correction to the particle propagator is the same 
as in the full theory [r| , and the ultrasoft correction van- 
ishes, so the wavefunction corrections are the same as in 
the full theory. The fermion graph, Fig. [5J gives 



Our normalization convention is such that the on-shell 
matrix element gets a contribution —5z/2 for each exter- 
nal particle. The wavefunction corrections for the vari- 
ous cases we need are tabulated in TableQl] In the table, 
we have distinguished between UV and IR divergences 
by the subscript on the 1/e terms. The scalar operators 
requires the scalar propagator correction, Fig. El which 
gives 
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FIG. 5: One-loop correction to the scalar propagator. 



The first term gives the wavefunction correction, and the 
second is the mass-shift of the scalar proportional to the 
gauge boson mass. The scalar mass-shift is canceled by 
the bare mass term in the scalar Lagrangian, which is 
adjusted to keep the physical scalar massless. This can- 
cellation is an example of fine-tuning required to have 
light scalars. 

The total on-shell amplitude /„ + I„ + I s — 2(Sz/2) is 



2 1 , 
- + -3- 



2U 



2L A/ L 



a/lq 



"kw - 3Lm 



5^ 
6 



(43) 



The total amplitude Eq. [[43]) is independent of S, n 
and r 2 introduced by the analytic regulator, and depends 
only on e of dimensional regularization. The cancellation 
of the v and S dependence is discussed in more detail 
in Appendix [XJ In evaluating Eq. (|43[) . we have used 



,+ - 



v 2 jp 2 and Q 2 = p{P 2 ■ The 1/e 



and 1/e 2 poles are ultraviolet divergences, and are can- 
celed by the renormalization counterterms in the effec- 
tive theory. The IR divergences in the EFT are regu- 
lated by the gauge boson mass, so the 1/e divergences in 
Eq. (|43|) are UV divergences. The 1/e term multiplied 
by —2 gives the SCET anomalous dimension listed in 
Table HI and is a non-trivial check on the analytic regula- 
tor computation. The SCET anomalous dimension was 
computed in Ref. [25[ using an off-shell regulator, and 
the analytic regulator gives the same result. While the 
total anomalous dimension is the same, the contribution 
of individual diagrams to the anomalous dimension de- 
pends on the regulator. For example, in Ref. (25|, the 
ultrasoft graph had a 1/e divergence which contributed 
to the anomalous dimension, whereas the ultrasoft graph 
vanishes on-shell when evaluated using the analytic reg- 
ulator method. Contributions can be moved between the 
collincar and ultrasoft diagrams, depending on the choice 
of regulator. 

The EFT below the matching scale /j, ~ M is SCET 
without gauge particles; thus there are no one loop di- 
agrams to consider in the theory below M. The fi- 
nite part of Eq. (|4"3"|) gives the multiplicative match- 
ing coefficient exp_D(/i) at the low scale /i of order M. 
The coefficient of O in the effective theory after inte- 
grating out the gauge bosons is given by c(/x — + ) = 



[exp -D(/i)] c(/i + + j. The coefficient D(n) has the usual 
expansion D = D^'a + D^ 2 'a 2 + . ■ ., and at one-loop 
order is 



-L 2 



2L M L Q - 3L 



m -r 77 



5^ 
6 



(44) 



for the fermion vector current. The other cases are com- 
puted similarly, and are given in Table |TJ The matching 
at M is independent of the vertex structure, and depends 
only on whether the particles are fermions or scalars. The 
-00 matching is the average of the results for two fermions 
and two scalars. This is a new feature of the effective the- 
ory, which follows because the graphs factorizc into con- 
tributions from the individual particles. The matching 
at Q from the full theory does not have this property. 

Note that D(/i) is a function of both I_m and Lq, and 
is linear in Lq . The matching condition depends on both 
scales Q and M. The dependence of the matching on the 
high scale Q is a new feature of SCET with massive gauge 
bosons, and has not occurred in previous computations 
in SCET, or other EFTs. We noted earlier that SCET 
graphs know about the scale Q through the labels n ■ p 2 
and n-pi since Q 2 = n-p 2 n-pi- Nevertheless, in previous 
computations such as DIS, the matching condition at the 
jet scale M 2 <C Q 2 depended on log M 2 /fj, 2 , and there 
was no log Q 2 /fi 2 dependence. It is easy to see why there 
must in general be Lq terms in the matching condition 
in our case. If D is the matching condition at [i, and "fh,i 
are the anomalous dimension in the theories above and 
below [i, 



d/i 



(45) 



In our example, "fi = since the theory below M is a free 
theory. Since jh bas the form Eq. (|3ip with a Lq term, 
such terms must also be present in D. Let us contrast 
this with DIS. For moments Mjv of the deep inelastic 
scattering structure function with N 3> 1, the jet scale 
is M 2 = Q 2 /N. The theory above the jet scale has an 
anomalous dimension 7^ which depends on Lq, and the 
theory below the jet scale has Altarelli-Parisi evolution 
with anomalous dimension 7/ which depends on log TV. 
The two anomalous dimension are related in such a way 
that 7; — 7h oc logQ 2 /7V, the logarithm of the jet scale. 
The matching D also depends only on the jet scale Q 2 /N, 
and there are no large logarithms in D if \i is chosen to 
be of order the jet scale [25| . 

The Lq term in Eq. ([4"4"|) is multiplied by La/, and so 
there is no Lq term in D if the matching scale is chosen 
to be exactly equal to M. This is accidental, and does 
not happen at higher orders. One can show explicitly 
that at two-loops, there is a non-zero Lq contribution to 
D even if /j, = M. In the standard model, it is convenient 
to integrate the weak gauge bosons out at a single scale 
fi = Mz, and one has one-loop terms in the matching 
condition of the form (logQ 2 /M|)(logM^/M|). 
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Renormalization group improved perturbation theory 
is used to sum logarithms in an EFT. This would not 
be possible if there were arbitrary powers of Lq in the 
matching condition. We will prove in Sec. EJ that to all 
orders in perturbation theory, the matching condition D 
is linear in Lq. Thus renormalization group summation 
can be used to obtain all logarithms except the first, so 
that in the Sudakov problem at order a n , the 2n — 1 
terms a n L 2n , . . . ,a n L 2 can be obtained by renormaliza- 
tion group evolution, but not the single log term a n L, 
which gets a matching contribution from D. The general 
form for D(fi) is 

D(n) = D (a(n),L M ) +D 1 (a(fx), L M )L Q , (46) 

which defines D ^. At one-loop, Eq. (j4"4")) gives 



D, 



(i) 



C F 



-L 2 m-3L 



9 5tt 2 



= 2C F l M . (47) 
Choosing [i — M gives the matching coefficient 

Q 2 



D(p = M) = D (a(M),0)+Z>i(a(M),0)log 



M 2 



(48) 



In our example, Di(a(M),0) — 0, so there is no 
logQ 2 /M 2 term in the one- loop matching coefficient. 
One expects that Di(a(M),0) ^ at higher orders, so 
there can be a single large logarithm in the matching 
coefficient. 

The final step in the computation is to compute the on- 
shell matrix element of O in the theory below M . Since 
the gauge bosons have been integrated out, this theory 
is a free theory, and the matrix element is trivial, being 
given by its free field value. The Sudakov form-factor 
is defined as the ratio of the scattering amplitude to its 
value in the free theory, so the low-energy matrix element 
contribution to the Sudakov form-factor is unity. 

The contributions to the Sudkakov form-factor are: 

1. The coefficient c(/i) in the full theory just above the 
matching scale = Q, which is chosen to be unity. 

2. The multiplicative matching coefficient exjpC(fi) 
for the matching between the full theory and SCET 
at the scale /x = Q. 

3. The integral of the SCET anomalous dimension be- 
tween \i = Q and fi = M . 

4. The multiplicative matching coefficient expl?(/i) 
for the matching at the scale [i — Q between SCET, 
and SCET with the gauge bosons integrated out. 

5. The low-energy matrix element, which gives unity, 
using the conventional normalization for the form- 
factor. 



Combining these contributions gives Eq. (|23f for the Su- 
dakov form-factor given earlier. The terms are repre- 
sented schematically as: 



C 7i 

Q — ► 



D 

M 



(49) 



The expression Eq. (f2"3"| for the Sudakov form-factor, 
with the one-loop coefficients given in Table Q] can be 
compared with known fixed order results in the case of 
the fermion vector current by expanding this in a power 
series expansion in a(M) as shown in Eq. (f2"5| . The result 
correctly reproduces the known cd_, a 2 L 4 and a 2 l_ 3 terms 
(L = logQ 2 /A/ 2 ). 

Comparison with the two- loop results of Ref. [7], [l5| al- 
lows us to extract values for the two-loop cusp anomalous 
dimension, 



+ — G F T F n s . 

The non-log part of the anomalous dimension is 
B( 2) = (4tt 2 - 3 - 48C(3)) C% 

-^ + 52 C (3)W. 



(50) 



- w + Y )C F T F n s 



(51) 



The log part of the matching at M for equal Higgs and 
gauge boson masses is 



(2) 



112 
~2T 



— 7T 
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C F T F nf 



782 
~27 



> 



c f 



where the Clausen function is 



Cl 2 (x) 



•s—^ smnx 



l 



n~ 



(52) 



(53) 



The anomalous dimension for Q > fi > M is indepen- 
dent of infrared physics, such as spontaneous symmetry 
breaking and the gauge boson mass, and so can be writ- 
ten in terms of group invariants such as C F and Ca ■ The 
expressions for and hold in a gauge theory with 
fermion and scalar fields in arbitrary representations. 

(2) 

The matching D\ depends on the gauge boson 
masses, and is only valid in a SU{2) gauge theory with 
scalars in the fundamental representation. The expres- 
sion Eq. (|52p has a C F T F n F term from fermion loop cor- 
rections to the gauge boson propagator, and a C F term. 
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FIG. 6: A graph whose group-theoretic factor cannot be writ- 
ten in terms of invariants such as CV and Ca- 



The Cp term arises from scalar loop corrections to the 
gauge boson propagator, as well as graphs such as Fig. [6] 
which arise due to spontaneous symmetry breaking. The 
group theory invariant for Fig. [6] depends on the pattern 
of symmetry breaking, and cannot be written in terms 
of SU(2) x U(l) invariants. Jantzen and Smirnov [l5| 
have therefore explicitly used the group theory factors 
for a broken SU{2) theory in evaluating these contri- 
butions, and we follow their convention here. Further- 
more, Ref. [l5[ computed the two-loop graphs only for 
Mh = Mw , and so Eq. (|5"2"j) is only valid for equal Higgs 
and gauge boson masses. 

The three-loop cusp anomalous dimension is known in 
a theory without scalar fields [Hf 



245 268 



44 



22 



-^T + ^ - irC(3) - t^tt 4 C 



27 



45 



836 80 2 112 „„A r T 



110 



- 32C(3) C F T F nf 



32 
27 



(T F n f f 



(54) 



so A^ 3 ) = 2r' 3 ' is known, neglecting scalar contributions. 
These missing scalar contributions are expected to make 
small corrections to T^ 3 ). The scalar term contributes 7% 
to the two-loop cusp anomalous dimension . 

Our one-loop computation combined with the known 
two-loop cusp anomalous dimension sums the LL and 
NLL series for the Sudakov form factor. The NNLL sc- 
ries requires A^ which is known excluding Higgs con- 
tributions, which is known (Eq. (|5Tj) L which 

is known for Mh = M, and which are known 

(Table |T|. For electroweak applications, the LL and NLL 
are more than adequate for precision studies. 



VI. PROOF THAT D IS LINEAR IN Lq 

The general functional form of the n-collinear graphs 
is expF(a(/x), L M , L 2 , L_), where L M = logM 2 /^ 2 , 

L 2 = logzyf/^ 2 and I = logu^/p^. Using = 

U i/Pi an( l Q 2 = P1P21 tins can be rewritten as 
exp F(a(n), Lm, L 2 , Li - Lq), where L Q = logQ 2 /^ 2 . 
Similarly, the n-collinear graphs have the functional form 
exp G(a(/i), l_Mj Li, L 2 — Lq). The sum of all the collinear 



graphs is the product exp(F + G), because the n and n- 
collincar graphs factor. These graphs give the matching 
coefficient expD, since the ultrasoft graphs vanish on- 
shell, so that D has the additive form 



D{a{fx), L M , Lq) 



F(a(fi), L M , L 2 , Li - Lq) 
+G(a(n), L M , Li, L 2 - Lq 



(55) 



The Li, L 2 dependence cancels, since D is independent of 

Equation (|55[) implies that -D(a(/i), Lm , Lq) is linear 
in Lq. The proof is as follows: Differentiating Eq. ([55]) 
with respect to Li and Lq, with respect to L 2 and Lq, 
and with respect to Li and L 2 give 



= d 1 d 1 F + d 2 d l G 1 
= d 1 d 2 F + d 2 d 2 G, 
= d 2 diF + d 2 d 1 G, 



where <9i. 2 is the derivative with respect to Li 2 . 
second derivative of D with respect to Lq is 



(56) 
The 



3D 
*5 



= dxdtF + d^G 

= -d 2 d 1 G-d 1 d 2 F 
= -d 1 d 2 {F + G) 
= -d x d 2 D = Q, 



(57) 



using Eq. (|56p and the commutation of partial deriva- 
tives, 8382 = 8283. Thus D can be at most linear in 
Lq. Equation (|55|) is D before the addition of renor- 
malization counterterms. The finite part shows that the 
matching correction is linear in Lq, justifying the form 
Eq. (|4"6"]l used earlier. The infinite part shows that the 
SCET anomalous dimension is linear i n Lq , and so gives 
another proof of this known result [H, [29| |. 



VII. CONSISTENCY CONDITIONS 

There are consistency conditions on matching coeffi- 
cients and anomalous dimensions which follow from the 
structure of the effective theory. Consider the match- 
ing of an operator between a high energy theory and a 
low energy theory at some scale y,. The operator coeffi- 
cients are Ch,i{/J>), with anomalous dimensions 7/i,z(/x) in 
the two theories, iidchj/d/i = jh,lCh,l- Assume that there 
is a multiplicative matching coefficient X(n) between the 
two theories, so that c;(/i) = X(/x)c/j(ju). The matching 
scale [i is arbitrary, so one gets the constraint 

At^-logQ = \i— logX + M-^-logc/j, (58) 
d// d/i d/z 



2 There are also wavefunction contributions to D. These are inde- 
pendent of Q. 
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which gives the relation 



Using Eq. (|46|) , and equating powers of Lq gives 



li - lh 



H—\ogX, 

a/j, 



(59) 



between the matching coefficient and the anomalous di- 
mensions in the two theories. In the Sudakov problem, 
applying Eq. (|59| to the matching between the full theory 
and SCET gives 



7' 



IF = V-T-C , 

d/_t 



(60) 



and applying it to the matching when the gauge boson 
is integrated out gives 



0-7 



(61) 



since the theory below M is a free theory and so has zero 
anomalous dimension. 

The anomalous dimension in the full theory and SCET 
have the form 7f(a(/z)) and 7 = — A(o(/x))Lq +B(a(/x)), 
respectively, where we have used the result that 7 is linear 
in Lq to all orders [25j |. The matching coefficient C at 
the high scale Q is independent of the low-energy scales 
such as M, and has the form C(a(/j), Lq). The matching 
coefficient D at [i ~ M can depend on both Q and M, 
since the SCET field labels depend on Q. As a result, 
D has the form D(a(fi), Lm, Lq). Any dependence on 
Q 2 /M 2 can be converted into dependence on Lm and Lq 
using Q 2 /M 2 = exp(L Q - Lm)- 



Equation (|60|) gives 



A(o(a*))L -B(o(/*))+7^(o(a*)) = 2 



Writing C as an expansion in Lq 



«9C 
5Lq~ 



(62) 



C(a(p),L Q ) = 53c n (o(/i))L5, (63) 



S(a( M )) = ^/? a (a)- 27^-2-^ 



9a 

A(a(/x)) = ^ a (a)-2|^. 
Writing D,;, i — 0, 1 as an expansion in Lm 



J <9L M ' 



(66) 



A(a(/i),L M ) = ^A,„(a(M))LM, (67) 



n=0 



gives the consistency conditions 
dD , 



- B = 



2nD 0>n = 



da 
dDn n 



P a - 2Di, - 2D ,i , 
-P a - 2Di, n - X , n > 2 . 



2nD 1:n (a) 



dD Un _ 



da 



-0a, n>2, 



(68) 



which determine -Di in , 72 > in terms of -E>i,0i A and 
B, and are satisfied by the one- loop values in Table |U 
The matching coefficient at fx = M is Dq o(a(M)) + 
D h0 (a(M))\ogQ 2 /M 2 . 



VIII. MASSIVE PARTICLES 

The calculations so far have been performed for ex- 
ternal particles with masses much smaller than the 
gauge boson mass. In this section, we extend the results 
to massive external particles. For fermions, we will as- 
sume that the theory is vectorlike, so that the fermion 
mass arises from a gauge invariant mass term — mipi/}. 
The standard model, a chiral gauge theory, where masses 
arise from Higgs couplings due to spontaneous symmetry 
breaking, is discussed in Sec. |Xj 



gives the consistency conditions 



dC 



lF {a)-B{a) = 2C 1 (a)-—^p a {a), 

oa 

df 1 

A(a) = AC 2 (a)-^-p a (a), 
dC 

2nC n (a) = —£±0 a (a), n>3, (64) 
oa 

which determine C n , n > in terms of Co, A, B and 7^, 
and are satisfied by the one- loop values in Table U The 
matching coefficient at n = Q is Co(a(Q)). 

Equation (fBTj) applied to the matching at M, gives 

A{a{p))L Q ~B{a^)) = ^/3 (a) - 2^- - 2^- , 

oa dig d\- M 

(65) 



A. Q > mi > M > mi 

Consider first the case where one particle has mass 777,2, 
with Q 3> 3> M , and the other particle has mass mi 
much smaller than M. For definiteness, the outgoing 
particle is taken to be heavier than M, and the incoming 
one lighter than M, but the results are symmetric under 
1 <-> 2. The Sudakov form-factor can be computed using 
a sequence of effective field theories. 

One first matches from the full theory to SCET with a 
massive particle at /1 ~ Q and uses the same set of oper- 
ators listed in Eq. (|30[) except £, n ,p 2 is now a n-collinear 
SCET field with mass m 2 as in Ref. [34]. This match- 
ing is independent of scales much smaller than Q, such 
as mi, m2 and M, and thus remains the same as in Ta- 
ble |U The second step is to run the operator in the 
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Pi 





(a) 



(b) 



FIG. 7: SCET graphs for the matrix element of O. The dotted 
lines are SCET propagators, and represent either fermions or 
scalars. The double lines are HQET propagators. 



--log 2 -^21og-f + - + 4 



(69) 



The wavefunction correction for a massive fermion is 

"12 ml 

+ — +4-31og— } 

_euv £ir M 



5z 



clCf 



(70) 



Combining Eq. (|69|) and Eq. ([70)) gives the multiplicative 
matching condition exp R,R = a + . . . , at [i of order 
m 2 , 



R 



clCf 



2 m2 



1. 



12 



(71) 



effective theory from Q to m%. The anomalous dimen- 
sion 7 is also independent of low mass scales and again 
gives the same result as the massless case. At the scale 
?7i2, one switches from SCET to a new effective theory 
in which the massive particle is described by a heavy 
quark effective theory (HQET) field h V2 , with a veloc- 
ity tj 2j with v\ = 1 [35|. The other (massless) particle 
is still described by a SCET field. The fermion vector 
current, for example, is now given by ^i, 2 7 M Wn£n,pi) in- 
stead of Eq. (|5D|) . and similarly for the other operators. 
The HQET field, h V2 , does not transform under collinear 
gauge transformations; therefore, there is no factor anal- 
ogous to the Wfj, Wilson line that goes along with the £ H 
field. h v still couples to ultrasoft gluons. One can make 
an additional field redefinition which eliminates the ul- 
trasoft gluon coupling to h v and introduces a Wilson line 
in the v 2 direction [36J. Both methods give the same 
on-shell matrix elements. 

The matching condition at 77i 2 is given by the differ- 
ence between the vertex graphs in Fig. [3] and Fig. and 
the wavefunction graphs, evaluated at /_t = 7772, i.e. be- 
tween graphs where £n )P2 and h V2 are used for particle 2. 
Note that there are three vertex graphs in the theory 
above 7772, and only two graphs in the theory below 7712, 
because there is no collinear Wilson line associated with 
h V2 . The graphs in the theory above 777,2 are evaluated 
with the gauge boson mass set to zero (since 7772 ^> M) 



at the on-shell point p 2 



The graphs in the ef- 



fective theory are evaluated with M — + at the onshcll 
point k% ■ v% = 0, where k 2 is the residual momentum of 
particle 2. Graphs Fig. [3Jd and Fig. [7^ are equal, since 
the n collinear interactions do not depend on whether the 
other field at the vertex is h v or ^ n ,p 2 ■ They cancel in the 
matching computation. The ultrasoft graphs Fig. [5J; and 
Fig. [TJd, the £n, Pl wavefunction graph, and the HQET 
wavefunction graph all vanish on-shell, so the matching 
computation is given by Fig. [3^, and the on-shell wave- 
function graph for £„ jP2 . The vertex graph Fig. [3Ji does 
not need an analytic regulator, and gives (for O = t/)7 m 7/>) 



1 



l/„ 777 

2 - log 7T 2 



where L m2 = log ml/ ' ■ The other cases are evaluated in 
a similar manner, and the results are summarized in Ta- 
blc lIIII The wavefunction correction for a massive scalar, 
which is needed for the last three rows, is 



Sz = aCp 



2 



2 



(72) 



The remaining steps are the evaluation of the anoma- 
lous dimension in the region M < \i < 7772, and the 
matching condition expS" at the scale M. These can 
be computed by evaluating the graphs in Fig. [7] and the 
wavefunction corrections, at the on-shell point &2 = 0, 
p\ = 0, and keeping the gauge boson mass non-zero. The 
graphs need to be regulated using an analytic regulator. 
The £n,pi is regulated as in Eq. (|33| for the collinear 
propagator. The Wilson line in Fig.[7K is regulated using 
Eq. (p4]) . as is the massless particle propagator in Fig.[7t>. 
The new feature is the HQET propagator for particle 2, 
which is regulated using 



1 



{-V2h) S2 
il+<5 2 



k 2 ■ v 2 [k 2 ■ v 2j 
Taking the HQET limit of the particle 2 propagator 



(73) 



2l l+fc 



[27772 k 2 ■ V 2 



l+<5 2 



(74) 



[(m 2 ?J2 + k 2 ) 2 — 777 

and comparing Eq. (fT3"|) with Eq. 
t/|/(2t772). Figure [7^, is given by the n-collinear graph 



we see that v 2 h 



evaluated earlier, Eq. 



2v H 



Figure [TJd is 

a, / -\ 6 



^ T(6 2 /2-S 1 /2)T(e + 6^2 + 6^2) 

The 77-collinear wavefunction graph is Eq. (I40|) . and the 
HQET propagator correction, Fig. [5] is 



clCf 



2 M 2 
2MTT+(k 2 -v 2 ) ( — + 21og — 
1 e /i z 



(76) 
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TABLE III: One-loop results for Q > 7712 > M > m 1 . R is the matching coefficient at fi ~ 7712, 72 is the anomalous dimension 
between m2 and M, and 5 is the matching coefficient at /i ~ M. The results only depend on whether the light particle is a 
fermion or scalar. 



k+p 




FIG. 8: HQET propagator Correction 



Eq. (|76[) gives a contribution to the heavy quark residual 
mass term, 5m — —laCpM-K^ and a wavcfunction con- 
tribution listed in Table HT1 The shift in the heavy quark 
mass is non-analytic in the gauge boson mass-squared. 
Such non-analytic contributions occur in mass correc- 
tions to particles with k ■ v propagators [13, . They 
arise from loop integrals which diverge as an odd power 
of k. Such integrals are finite, but non-analytic, in di- 
mensional regularization. 

Adding Eqs. (|3"5j) and Eq. (|75|) . expanding in 5i, 
and subtracting the wavefunction corrections due to 
the heavy fermion, Eq. (|76[) . and the collinear fermion, 
Eq. dU]), gives 



aCf 



1 



1 



- + - - + 21og 



fjLm,2\ 



9 5tt 2 
4 ~12 



2 M 



5 ^ 
2 l0g M^ 

/_tm 2 



- - — l0g Z -^r + 2 log At lOE 



Q 2 



(77) 



To obtain Eq. ([77]) we have used v 2 h = v \l^ m i)i v \ = 
vf/Pi, v% — ^2/^2" > anc ^ Q 2 — m 2Pi ■ The 1/e coefficient 
multiplied by —2 gives the anomalous dimension 



72 



aC F (-5 - 2L r , 



4L 



Q) 



(78) 



and the finite part gives the matching correction 

5. 



,5' 



aCf 



-M 



9 5tt 2 
4 ~ "l2~ 



FIG. 9: Vertex Correction in the theory below 7711,2- 



"2 L m ~ l m (L„ i2 - 2Lq) 



(79) 



The other cases are computed similarly, and are given in 
Table HU 

The terms which contribute to the final result are sum- 
marized schematically as: 



C 71 R 72 5* 
Q — > m 2 — > M 



(80) 



B. 



»m 2 »mi» M 



The second case we consider is where both particles 
have mass between Q and M. For definiteness, we choose 
7712 > ™>i- The Sudakov form-factor can be computed 
using a sequence of matching and running steps. The 
matching expC at Q, the anomalous dimension 71 be- 
tween Q and 77J2, the matching expi? at 7772, and the 
anomalous dimension 72 between 7772 and mi are all in- 
dependent of the lower scales mi and M, and have the 
same values as in Tables U IIII1 

The new feature is the matching condition exp T at the 
lower particle mass mi. The graphs in the theory above 
mi are shown in Fig. [7j In the theory below mi, the 
SCET field for particle 1 is replaced by the HQET 
field h Vl . The fermion vector current, for example, is now 
given by h V2 j^h Vl instead of hy^^W^n^ ■ The vertex 
correction in the theory below mi is shown in Fig. [9] 
There is only one vertex graph instead of two, because 
there is no collinear Wilson line W associated with the 
HQET field h Vl . The matching condition is given by 
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TABLE IV: One-loop results for Q > m,2 > m\ > M. T is the matching at mi, 73 is the anomalous dimension between mi 
and M, and (7 is the matching at M. T only depends on whether the light particle is a scalar or a fermion 



computing the difference of graphs Fig. [7] and Fig. [9] on- 
shell, and setting all scales less than mi to zero. The only 
non-zero contribution is from Fig. [7^ and the fi-collinear 
wavefunction renormalization graph. These are the same 
graphs that contribute to the matching condition at m.2, 
so T is given by R with mi — > mi, and is tabulated in 
Table [TV] 

The remaining quantities needed are the anomalous di- 
mension 73 between m\ and M , and the matching condi- 
tion exp U at M . These can be computed from the graph 
in Fig. [9] evaluated on-shell, but now with the gauge bo- 
son mass M included. The graph gives 



a wr(w) 



where w = V2 ■ v\ and 



2 M 2 
h 2 log — 2" 



r(w) 



log (w + yt? 



(81) 



(82) 



is the factor which occurs in the velocity-dependent 
anomalous dimension in HQET 35|. Including the 
heavy quark wavefunction correction, Eq. (|76p . gives the 
anomalous dimensions and matching coefficient listed in 
Table HVl In the high energy limit Q 2 ~ 2m,im2W, 



wr(w) - log(2iu) . (83) 
The terms which contribute to the final result are: 



C 71 R 72 T 73 U 

Q — ► m,2 — > mi — > M 



(84) 



C. Q > m 2 = mi » M 

The third case we consider is where the two particles 
are degenerate, with m-2 = m\ = m, and Q » m » M. 
This can be computed using the results already derived. 
The matching exp C at Q and the running 71 between Q 



3 T depends on whether the particle being 
fermion or a scalar. 



integrated out is a 



and m is the same as in Table |TJ The matching at m is 
given by switching both particles from SCET to HQET 
simultaneously. This is just the sum of the matching 
coefficients at mi and m\ computed previously, so the 
matching condition is exp [R + T], with mi = mi = m, 
where R and T are given in Tables IIIII and HVi respec- 
tively. The anomalous dimension between m and M , and 
the matching at M are given by 73 and U in Table HVl 

The terms which contribute to the final result are: 



C 71 
Q — » 



R 

m 2 - 



T 73 
mi — > 



U 
M 



(85) 



D. Q > m 2 ~ M > mi 

It is also useful to derive results where the particle 
masses and gauge boson masses are not widely separated 
from each other. In this case, it is more important to 
include the full m/M dependence, rather than sum high 
order a log m/M terms, which are no longer very large. 
In the standard model, this situation arises for the W, 
Z and t, which have masses which are not sufficiently 
widely separated that electroweak logarithms need to be 
summed. If two (or more) particle masses are not widely 
separated, one can make a transition to a new EFT by in- 
tegrating out both particles at a common scale /x, rather 
than integrating them out sequentially. The results for 
the various cases are summarized in this and the follow- 
ing subsections. The difference from previous results is 
that one has to include all the relevant masses in the par- 
ticle propagators. For example, for the case studied in 
this subsection, m2 ~ M, one also includes m.2 in the 
denominator of Eq. (|32p . The integrals now depend on 
a dimensionless parameter, the ratio of particle masses 
m\jM 2 . 

If Q ^> m2 ~ M mi, then the matching at Q and 
the running between Q and 7712 ~ M, remain unchanged, 
and are given by C and 71 in Table HI At the scale fi of 
order m-2 ~ M, one integrates out the gauge boson, and 
switches to a theory in which particle 2 is described by 
a HQET field. In this matching, the n-collinear graph 
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Eq. (|32]1 with the analytic regulator is now 



I n = ~ig 2 ^C F c{^i) 



d d k 



1 



{2n) d k 2 - M 2 2 
Si 



2 [(pa - fc)2 - ml] 1+52 hrl-fc] 



and is evaluated at the on-shell point p 2 



l+5i 



(86) 



The 



collinear and ultrasoft integrals remain unchanged. We 
saw in Sec. IVIII Al that the S and v dependence can- 
celed between all the diagrams for the massless case. The 
cancellation must still hold when /„ is evaluated with 
m 2 7^ 0, so that aC F f F (z) = I„(m 2 ) - I n (m 2 = 0), 
z = m 2 /M 2 has a finite limit independent of the analytic 
regulator, as can be verified by explicit computation. The 
result for f F (z) is given in Eq. (|B1[) . The wavefunction 
renormalization is also modified, and the shift hp(z) is 
given in Eq. (|B2|) . The matching condition D is given by 
D{m 2 ) = D(m 2 = 0)+aC F (/f(*) - h F (z)/2) if the par- 
ticle integrated out is a heavy fermion, where the mass- 
less value D(m 2 — 0) is given in Tablc|TJ Similarly, if the 
particle integrated out is a heavy scalar, the matching D 
is D(m 2 ) = D(m 2 = 0) + aC F (fs(z) ~ hs{z )/2) where 
the scalar functions are given in Eqs. (|B4IB5[) . 

The theory below m 2 ~ M is a theory in which parti- 
cle 2 is described using a HQET field, and particle 1 by 
an SCET, with the massive gauge boson integrated out. 
In our toy example, this is a free theory. 

Schematically, the terms are (z = m\/M 2 ): 



C 7l D + aC F (f(z)-h(z)/2) 
Q — ► m 2 ~M 



(87) 



the same as for m, = 0. The n and n collinear graphs at 
the scale mi ~ m 2 are independent of each other, so the 
matching is given by R(m 2 ) +T(mi) given in Tables HTT1 
IIVI Below mi ~ m 2 , the computation reduces to that in 
Sec. IVIII CI 

Schematically, the terms are: 



C 7i 



R 

mi ' 



T 73 
m 2 — ► 



G. Q>rm 



TO2 



u 

M 
M 



(89) 



The evolution to /i ~ mi ~ m 2 ~ M is the same as for 
the massless case. The matching at /i involves massive 
collinear propagators, each of which modifies the massless 
matching condition, so the matching is given by 

D{m 1 ,m 2 ) = D{m l = m2 = $) + aC F {h{z 2 )-h 2 {z 2 )/2) 
+aC F (f 1 (z 1 )-h 1 (z 1 )/2) (90) 

where Zi — mf/M 2 , and f\ y2 , hi >2 are chosen to be 
fF.s and hp,s depending on whether the correspond- 
ing particle is a fermion or scalar. The massless value 
D(m\ = ra 2 = 0) is given in Table |TJ 

Schematically, the terms are [z^ — mf/M 2 ): 

C 7i D + aC F (f(z 1 )-h(z 1 )/2) + aC F (f(z 2 )-h(z 2 )/2) 
Q — > m\ ~ m 2 ~ M 



(91) 



IX. SCALAR CORRECTIONS 



E. Q > m 2 > M ~ mi 

If Q m 2 3> M ~ mi, the matching and running 
down to M ~ mi remains the same as in Sec. IVIII Al 
In the theory at M <~ mi, particle 2 is described by an 
HQET field, and particle 1 by a n collinear field. The 
matching condition at M ~ mi would be given by S if 
mi — > 0. By the same arguments as above, the effect 
of mi is to modify the fi-collinear integral by a finite 
amount, so the matching is now S + aC F (f(z) — h(z)/2), 
with z — ml/M 2 , where /, h are the fermion or scalar 
values Eqs. (|B1IB2|) or Eqs. (|B4IB5|) . depending on the 
type of particle 1. 

Schematically, the terms are [z = m\/M 2 ): 

C 7i R 72 S + aC F (f(z)-h(z)/2) 
Q — ► m 2 — > mi ~ M 

F. Q > mi ~ m 2 > M 

The situation is similar to Q 3> mi = m 2 3> M consid- 
ered in Sec. IVIII CI The evolution down to mi ~ m 2 is 



In this section, we compute the scalar exchange cor- 
rections to the Sudakov form-factor. The graphs are the 
same as those for gauge exchange, with the gauge boson 
replaced by the scalar %, with mass M x . As for gauge 
bosons, one needs to include both collinear and ultrasoft 
fields for \ to represent collinear and ultrasoft x parti- 
cles. In the gauge boson results we removed an overall 
factor of a = g 2 / (16ir 2 ) . In the %-exchange graphs, we re- 
move a factor of hih 2 /(16n 2 ), where h% >2 are the Yukawa 
couplings at the two vertices. Note that the coupling of 
\ to scalars h^j has dimensions of mass, so the factor 
removed for operators such as <jy<j> is dimensionful. Un- 
like for gauge interactions, the wavefunction and vertex 
corrections can have different coupling constants. 

Collinear gauge bosons and matter fields in the op- 
erator O occur in the combination W\^ n ^ p or W£<S> niP . 
They arise from full theory graphs in which the gauge 
fields couple to the other particle, which moves in the 
fi-dircction. The intermediate propagators are off-shell 
by order Q 2 , and can be shrunk to a point, as shown 
in Fig. [101 Single gauge boson emission gives the ver- 
tex gn^/^—n ■ k). At higher orders, multiple gauge bo- 
son emission from the n particle line combined with the 
non-Abelian multi-gluon interaction give the Wilson line 
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FIG. 10: Graphs in which n-collinear gauge bosons couple to 
the n-collinear line generate the Wilson line W„. 



operator W n . Multi-gluon emission is related to single 
gluon emission by gauge invariance, and this relation 
holds even in the presence of loop-corrections. The Wil- 
son line structure of the vertex W n is required by collinear 
gauge invariance [HI]. 

One has a similar construction for multiple n-collinear 
X fields emitted from the n particle. At the level of single 
X emission, the x vertex is h^ t ^/(—n ■ k) instead of the 
gauge vertex gn^ / (— n ■ k). This is all we require for our 
computation. It would be interesting to work out the 
structure of the scalar vertex at higher orders, including 
radiative corrections. Multi-scalar emission is not related 
to the single scalar vertex by gauge invariance. 

Most of the scalar corrections vanish. In SCET, the 
fcrmion Yukawa vertex vanishes, because Eq. ([3|) implies 
that 



(92) 



using i/n/i = n 2 = 0. The tri-scalar couplings x</>^</> have 
dimensions of mass, and x exchange corrections to the 
scalar operators are suppressed by powers of h^/Q, which 
is subleading in the EFT power counting given our as- 
sumption that does not grow with Q. The easiest way 
to see this is to use the rescaled fields 4> n ,p, which have 
a propagator of the same form as that for fermions. The 
Yukawa coupling becomes 



Il4 

n ■ p 



■ X': 



n,pYn,p j 



(93) 



which is order l/Q since n • p is order Q, and gives an 
explicit l/Q suppression to the graph at each tri-scalar 
vertex. 

One interesting point is the decoupling of scalars be- 
low mi and ni2, so that particles 1 and 2 can be treated 
as HQET fields. If one directly matches from the full 
theory onto HQET, then h^x^P^ ~^ h,pxhvh v , which 
is non-zero. If instead, one first goes through SCET, 



then h^xiPi' ~> hil>x€n,p€n 



0, so we have an apparent 



contradiction. However, the two results are not in dis- 
agreement. The scalar HQET vertex graph (Fig. [5] with 
the gauge boson replaced by a scalar) is equal to —1/w 
times the corresponding gauge graphs, rescaled by the 
ratio of the Yukawa couplings to the gauge couplings. 
In the Sudakov limit, w ~ Q 2 /( m i m 2)j an d \/w is a 
power suppression which can be neglected, so both ways 
of matching agree, since power corrections are neglected. 
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TABLE V: One- loop scalar exchange contribution to on-shell 
wavefunction renormalization. The exchanged scalar mass is 
M x , L x = log M x /fj, 2 , and the particle (fermion or scalar) 
mass is m. /if,s are given in Appendix[B] An overall factor 
of hih 2 /(l6TT 2 ) is omitted. 



The only scalar graphs which remain are the match- 
ing at Q, which are full theory graphs, and scalar con- 
tributions to wavefunction renormalization in the effec- 
tive theories. The wavefunction contributions are sum- 
marized in Table \V\ where L x = log M 2 /[i 2 . The full 
theory wavefunction renormalization vanishes for both 
fermions and scalars, so the entire matching correction 
at Q given in Eq. (|94|) arises from the vertex correction. 
The EFT matching and running can be computed from 
the scalar wavefunction graphs in Table fVl The EFT ma- 
trix elements are given by taking —1/2 times the entries 
in the table for each particle, multiplying by h 2 /(16ir 2 ), 
and then adding the contributions from the two parti- 
cles. The finite part gives the matching correction, and 
(—2) times the coefficient of the 1/ctjv term gives the 
anomalous dimension. Since only wavefunction graphs 
contribute, there are no Lq terms which can only arise 
from vertex graphs. 

The computation of scalar contributions to the anoma- 
lous dimension and matching for the various cases con- 
sidered in Sec. I Villi parallels the gauge boson discussion. 
The matching coefficients for mi 2 =/= are given by a 
formula analogous to Eq. (|90p. with the gauge boson 
functions /i?s and h,F,s replaced by the corresponding 
X-exchange functions fp,s ~ * 0, since there are no vertex 
corrections in the effective theory, and fiF,s ~ * hp,s- 

We have divided the scalar exchange contributions into 
vertex and wavefunction pieces, rather than giving the 
total contribution as in the gauge case. The reason is 
that the standard model is a chiral gauge theory, and the 
Yukawa couplings connect one matter representation to 
another. Thus, vertex corrections can mix left-handed 
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currents with right-handed currents, whereas wavefunc- 
tion corrections do not mix different SU(2) x U(l) rep- 
resentations. Keeping the two contributions separate al- 
lows us to compute Higgs radiative corrections in the 
standard model using the the results given in this sec- 
tion. 

The matching corrections C at Q can be computed 
as for the gauge boson case. The matching corrections 
for external scalar particles due to scalar \ exchange 
are power suppressed, and vanish to leading order. For 
fermions, the x exchange corrections give 



iptp 



exp [a Y (-2 + L Q )] [^, P2 W n }[W^ n , Pl } 



exp 



(n"n" - n v n") [U P2 W n }[wtn, Pl . 



ay = 



16tt 2 



(94) 



which have to be combined with the gauge boson match- 
ing conditions in Table [T] The vertex graph contributes 
2ay, —ay and 0, respectively, to the anomalous dimen- 
sions of the three operators in the full theory. The wave- 
function graphs contribute an additional ay to all three 
operators. 



X. APPLICATION TO THE STANDARD 
MODEL 

The results we have obtained for the toy theory can 
now be used to compute results for the standard model. 
One has to be careful in using the correct coupling con- 
stants, since the standard model is a chiral gauge theory, 
and the toy model is vector-like. 

For the LHC, one is interested in processes such as di- 
jet production. The SCET operators at the high scale Q 
involve more than two SCET fields. E.g. in qq — > qq, the 
EFT operator has four fields, two for the incoming parti- 
cle and two for the outgoing ones. One can obtain results 
for more than two external particles by combining the 
two-particle results computed in this paper with the ap- 
propriate gauge theory factors such as Ca and C p ■ There 
are several interesting features of the analysis which are 
independent of the calculations presented in this paper, 
so we defer the discussion of experimentally relevant ex- 
amples to a subsequent publication [20l j . Here we show 
how our results can be used to compute the radiative cor- 
rections to quark production by a gauge invariant current 
Qd^PhQii where Qi is the quark doublet 4 for generation 



i = u,c, t, and to charged lepton production by Lj^PlL. 
We will do the computations for light quarks in Sec. IX A[ 
for leptons in Sec. IX B[ and for top quarks in Sec. IX CI 
All fcrmion masses other than the top quark mass are 
neglected. 



A. Light Quarks 



The first generation quark doublet is 



Qu 




Vudd + V us s + Vubb 



(95) 



using the mass eigenstate basis. At the scale Q » m q 
the coefficient of the operator in the full electroweak the- 
ory is assumed to be unity. For the first generation, all 
quark masses and Yukawa couplings can be neglected, 
and so the answer is given by combining the gauge boson 
contributions computed earlier. 

The operator in SCET at the scale Q is 



Qu^PlQu - <Q)&W n ]^P L [wU^], 



(96) 



where represents the left-handed electroweak u- 

quark doublet Eq. (|95|) in SCET, and we have suppressed 
gauge indices. The matching condition is 



log c(Q) = 



q s (Q) 4 o 2 (Q) 3 mjQ) 1 
4tt 3 4?r 4 4tt 36 



(97) 



using the third column of Table [T] with Lq = at the 
scale fi = Q. The gauge couplings have been multiplied 
by the corresponding Cf values: 4/3 for an SU(3) triplet, 
3/4 for an SU(2) doublet, and 1/36 for Y = 1/6. The 
electroweak couplings renormalized at fi = Mz are 



a 2 {M z ) 
ui{M z ) 



«em(Mz) 

sin 2 9w(M z ) ' 

a C m(M z ) 
cos 2 e w (M z ) ' 



(98) 



and their values at Q are obtained by the usual /3- 
functions of the standard model. 

The theory below Q is SCET with an SU(3) x SU(2) x 
C/(l) gauge symmetry. In this regime, the SCET current 
in Eq. (|96|) is multiplicatively renormalized with anoma- 
lous dimension (from the fourth column in Table HJ 



7(m) = 



o= s (m) 4 
4tt 3 



47T 4 



ai(jj) 1 



4tt 36 



[4Lq - 6] 



(99) 



4 Not to be confused with Q, the momentum transfer. 



The anomalous dimension 7 is used to run c down to 
a scale of order the gauge boson mass. One can inte- 
grate out the weak gauge bosons sequentially, by first 
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integrating out the Z boson at fx = Mz , followed by the 
W at fi = M w . This sums (a log 2 M w /M z ) n , n > 1 
terms, while neglecting a(Mw/Mz) n , n > power cor- 
rections. This is not a good choice to use for the standard 
model, since Myy/Mz is not very small, and summing 
powers of M\y/Mz is more important than summing 
a log 2 Mw/Mz terms. Instead, we integrate out the W 
and Z at a common scale, chosen to be \x = Mz- In this 
way, we match directly from an SU(3) x SU(2) x U(l) 
gauge theory onto a SU(3) x U(l) cm gauge theory of 
gluons and photons, and there are no complications of 
an intermediate stage of broken electroweak symmetry 



where the Z is integrated out, but not the W. 

At the scale fi — Mz, integrating out the W and Z 
bosons give a matching correction to the SCET operator, 

[^W n ]-fP L \Wi^] - a^[^l 2 W n ]^P L [Wk^l] 
+a id ' ) [&i\W n ]rPL[WkS 1 ] ■ (100) 

Since the electroweak symmetry is broken, the u and d! 
parts of the operator get different matching corrections. 
The corrections dj are obtained using the last column of 
Table ED 
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(101) 



The first term for log 

a (u,d') i s 

the Z contribution, the 
second term is the W contribution, and the coupling con- 
stants are renormalized at Mz- 

Below Mz, the operators in Eq. (jlOOp are multiplica- 
tively renormalized, with anomalous dimensions 



[4L Q - 6] , 

[4L Q - 6] , (102) 
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The final result for the operator at a low scale is 

+c ( ^[^W„]7 ,, i , L[W'M?J ) 
with 

' Mz d M 



(103) 



logc^(/x) = log c(Q) + 



+iog w+ r ^7 (u) (m), 



logc^V) = log C (Q) 



d/x 
A* 



t(m) 



loga (d,) + 



" *W),(104) 



where the various pieces are given in Eqs. (|97l99ll01ll02p . 
The EFT operator Eq. (|103|) can then be used to compute 
processes such as dijet production using SCET [39(. For 
jet production, the scale \i would be chosen to be of order 
the jet invariant mass, around 30 GeV for jets at the 
LHC. 



B. Leptons 

The computation for the radiative corrections to the 
lepton current Lj^P^L, where L is the lepton doublet 



L = 



(105) 



is similar to that for the quark doublet, and we summa- 
rize the final result. The full theory operator at the low 
scale n is 

L^P L L -> c^[^ P2 W n }^P L [Wt ( nL] 
+c {1 \^MYPl[WU^ p1 ] , (106) 



with the coefficients given by Eq. (|104[) with u — > v, 
d! — ► £, where the on the rhs of Eq. (|104[) for leptons are: 



log c(Q) = 



q 2 (Q) 3 ai(Q) l 
Air A At: A 
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(107) 
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C. Top Quarks 

In this subsection, we show how our results can be used 
to compute the radiative corrections to tt production by 
a gauge invariant vector current Q^PlQu where Q t is 
the left-handed quark doublet in the standard model, 



Qi 




V td d + V ts s + V tb b 



(109) 



and b' — V tc [d + V ts s + V t bb using the mass eigenstate 
basis. We will neglect all quark masses other than m t - 
This example illustrates how to use the fermion mass 
and Higgs exchange contributions computed in the toy 
example. 

The operator in SCET at the scale Q is 



- ci(QM ( n Q ^w n }rPL[wM 



(Qth 
Pi 1 



(110) 



where fW*) and represents the left-handed elec- 
troweak t-quark doublet Eq. (|109[) and the right-handed 
t-quark singlet £r in SCET and we have suppressed 
gauge indices. The tji terms arise from Higgs exchange 
graphs Fig. [TT] 

The matching condition is 



logci(Q) 
c 2 (Q) 



a s (Q) 4 t a 2 (Q) 3 

47T 3 47T 4 
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16tt 2 
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[111) 



using Tables U and Eq. (|94|) with Lq = at the scale 
fj, = Q. The gauge couplings have been multiplied by 
the corresponding Cp values, 4/3 for an SU(3) triplet, 
3/4 for an SU(2) doublet, and 1/36 for Y = 1/6. 




FIG. 11: Higgs correction which causes Qt7 M PhQt to mix 
with tRj^PRt. The index a is an SU(2) index, and is summed 
over. 



The top quark Yukawa coupling is normalized so that 
gt = \f2rnt/v, with v ~ 247 GeV. The Higgs exchange 
graph in the chiral standard model has been computed 
using the toy-model value Eq. (|94|) for a vector-like the- 
ory, combined with the result that a Yukawa vertex flips 
the fermion chirality. The factor of 2 in front of <? 2 /(167r 2 ) 
arises from summing over a closed SU(2) index loop, i.e. 
because the Higgs and Q are 5/7(2) doublets (the sum 
on a in Fig. [TTj) . The factor of 1/2 in square brackets is 
the coefficient of ay in the second line of Eq. (|94|) , with 
Lq — > at — Q. The Higgs exchange vertex correc- 
tion mixes the Ql operator with the tn operator. Higgs 
exchange corrections do not contribute to the diagonal 
coefficient c\ , since the first row of Table [V] shows that 
the full theory wavefunction renormalization has no finite 
part. 

The theory below Q is SCET with an SU(3) x SU{2) x 
U(l) gauge symmetry. In this regime, the two opera- 
tors in Eq. (| 1 1 0[) are multiplicatively renormalized with 
anomalous dimensions 
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The ti, wavefunction factor due to Higgs exchange does 
not have the factor of two from the SU(2) index sum- 
mation that is present for tn- The Higgs vertex graph, 
which causes c\ — ci mixing, is 1/ Q 2 suppressed. 

The anomalous dimension 7 is used to run down to 
a scale of order mt ■ At this scale there are several differ- 
ent methods one can use. As for massless quarks, one can 
integrate out the scales m t , Mw, Mz and Mh in various 
ways, e.g. one can integrate out each particle at a scale /i 
equal to its mass, or integrate out one or more particles 
simultaneously at some common value of /1. Integrating 
out the top quark leads to a complicated effective the- 
ory with dynamical W and Z bosons which is no longer 
SU(2) x U (1) invariant, since the V quark is in the theory 



but not t. Luckily, the best method for experimentally 
relevant computations is also the simplest to use: Since 
m t , Mw, Mz, and presumably Mh are not widely sep- 
arated, one can integrate them all out together. In this 
way, one goes directly from an SU(3) x SU(2) x U(l) 
invariant theory to a SU(3) x £7(l) em gauge theory, with 
broken SU (2) x !7(1) symmetry and no electroweak gauge 
bosons. This procedure keeps the entire mass dependence 
on the four mass scales. 

At the scale /j, = mt the f-quark SCET field is replaced 
by the heavy quark field t v , whereas the b' quark SCET 
field in the doublet fW*) remains an SCET field £ (b '\ 
The operator matching is 



1 - 

— &lt v2 t vl 



-Q3t V2 t Vl , (113) 



where the matching coefficients are denoted 01.2,3. Using 
the result of Sec. IVIII CI for the gauge boson exchange 
graphs, and Sec. [IX] for the Higgs exchange graphs gives 
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All running couplings are renormalized at fi = m t . The 
expressions are given by adding the contributions due to 
the Z (F g (Q,Mz,m t ) term), W {FJQ,M w ,m t ) term), 
gluon, 7, H (F h (M H ,m t ) term), h° (F h (M z ,m t ) term) 
and h + (Fh(Mw,mt) term), where h°, h + are the un- 
physical Higgs scalars present in R^—\ gauge. 



Below (i — mt, the t V2 t Vl operator has anomalous di- 



mension (from the third column of Table IIV[) 
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where w = V2 • V\ — 1 + Q 2 / (2m 2 ). 

The radiative corrections to the it operator can then be 
combined with known methods to obtain t-quark decay 
distributions [36[. The QCD corrections (the a s terms) 
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have already been included in the analysis of Ref. [36|. 
The new results in this paper are the additional elec- 
troweak radiative corrections, including Higgs effects. 



XI. NUMERICS AND CONCLUSIONS 

We have shown how SCET methods can be used to 
compute the radiative corrections to electroweak pro- 
cesses. We discussed the results for massless external 
particles given in Ref. 19] in more detail, and derived 
the result that there is at most a single power of Lq in 
the matching at M to all orders in perturbation theory. 
The existence of Lq terms in the matching at M is a new 
feature of SCET with massive gauge bosons. The re- 
sults of Ref. have been extended to include external 
particle masses and radiative Higgs exchange corrections 
proportional to the Yukawa couplings. 

Most of the paper used the vector-like SU{2) gauge 
theory. Section [X] explained in detail how the the re- 
sults for the vector-like theory could be used to compute 
radiative corrections in the standard model, which is a 
chiral gauge theory. In this paper, we have computed 
radiative corrections to operators with two external par- 
ticles. These can be used to compute the production 
rate for two external particles by a gauge invariant bilin- 
ear source. This could be applied to the decay rate of a 
(hypothetical) gauge singlet particle into two fermions. 
SCET methods can also be used to compute the radia- 
tive corrections to the dominant high-energy processes 
observable at the LHC, such as dijet production from 
quark-quark scattering. For a partonic process such as 
qq — ► qq, the EFT operator is a four-quark operator, 
and the radiative corrections can be obtained by using 
the results of this paper, summed over all pairs of parti- 
cles. The anomalous dimensions are about twice as big 
as the ones for the two-quark operators considered here. 
We postpone further discussion to a subsequent publica- 
tion Q. 

We conclude by giving some plots which show the typ- 
ical size of the radiative corrections. The LHC center 
of mass energy is y/s = 14 TeV. The partonic center of 
mass energy y/I is much lower, because a proton with 
energy E has partons with energy fraction x < 1, given 
by a parton distribution f(x) which vanishes as x — ► 1. 
The bulk of the dijet cross-section is at low s, but the 
LHC has sufficient luminosity to be able to observe the 
dijet cross-section up to y/1 of order several TeV. We 
plot the Sudakov form factor for electron production via 
Lj^PlL, it-quark production via QuI^PlQu and i-quark 
production via Qtl^PhQt, which are the results given in 
Sec. |Xj for Vs between 0.25 and 8 TeV. Figure [121 gives 
the results for Fe(Q 2 ) for the three cases, stopping the 
evolution at \i = Mz- In Fig. [T31 the EFT operators 
have been evolved all the way down to /i = 30 GeV, 
the typical invariant mass used to define a jet at the 
LHC. Figure [TJ] shows the electroweak contributions to 
the three form-factors as a percentage change relative to 
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FIG. 12: The Sudakov form-factor for it-quarks (solid black), 
t-quarks (dotted red) and electrons (dashed blue) at fj, = Mz 
for m H = 200 GeV. 
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FIG. 13: The Sudakov form-factor for it-quarks (solid black), 
t-quarks (dotted red) and electrons (dashed blue) at fi = 
30 GeV for m„ = 200 GeV. 



the form factor including only the QCD radiative cor- 
rections. We have used a Higgs mass of 200 GeV in the 
plots. Varying the Higgs mass between 150 GeV and 
500 GeV makes a difference of less than 0.5%. Fe 
is normalized to unity in the absence of radiative cor- 
rections. Radiative corrections for electrons are about 
7% at 2 TeV, and \i = Mz, increasing to about 8% for 
H = 30 GeV due to the QED running below M z . The 
corrections for quarks are much larger, 15% for the t- 
quark, and 30% for the it-quark at /i = Mz, increasing 
to 40% and 50%, respectively, at \i = 30 GeV. There is 
also a significant difference between the results for t- and 
it-quarks, arising from power corrections which depend 
on mt/Mw,z and Higgs corrections which depend on the 
Yukawa coupling g t . The bulk of the difference is due 
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FIG. 14: The electroweak contribution to the Sudakov form- 
factor (as a percentage change) for u-quarks (solid black), 
t-quarks (dotted red) and electrons (dashed blue) at p — 
30 GeV for m„ = 200 GeV. 



APPENDIX A: ANALYTIC REGULATOR 
DEPENDENCE OF COLLINEAR GRAPHS 

The n-collinear contribution Eq. (I37|) can be written 



where 
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to the power corrections, the gt terms are less than 1%. 5 
The corrections for the four-quark operators needed for 
realistic processes are bigger, with corrections even to 
color singlet processes being greater than 20%. The ra- 
diative corrections are large enough that resummation is 
necessary to get an accurate prediction for the partonic 
cross-sections. We have shown how one can perform the 
resummation using EFT methods. Previous computa- 
tions 0, 0, H 0, 0, 0, 1, 0, M, d d, El Q d have done 
the resummation using infrared evolution equations 
The EFT method allows one to include mass effects as 
well as Higgs corrections in a systematic way, which 
have not been included previously. It also handles the 
cross-over between the SU(3) x SU(2) x U(l) and the 
SU(3) x U(l) gauge theories above and below the weak 
interaction symmetry breaking scale, including the ef- 
fects of unequal W and Z masses. The infrared evolution 
method uses a conjectured form for this cross-over with 
equal W and Z masses. The extension of previous results 
to massive external particles is currently being studied by 
other groups using a method-of- regions analysis [401 ] . 
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The n-collinear contribution is given by Eq. (|A1[) with 
Tj -> -r), p -> -p, vi -> vi, Vi/pz -> v£ /pf. The 
dependence on n,2 and v\,2 is additional scheme depen- 
dence introduced by the analytic regulator. The sum 
of the n-collinear and n-collinear contributions simpli- 
fies greatly, and gives Eq. (|4*5|) on using v± = vf/pf, 



v. 



— v'i/p^- In particular, the n and p terms cancel 
between the two contributions. The relevant relations 
are 
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The identities Eq. (|A3|) do not make any assumptions 
about the values of v\ and v 2 , which need not be equal. 

It is convenient to use the special case of Eq. (|Alj) with 
p = n = to define the n, n-collinear contributions. This 
form is given by using the analytic regulator followed by 
the limit r 2 = — ri, ri — > oo, and gives 
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5 Note that the Higgs corrections do not have the Lq enhancement 
that is present for the gauge bosons. 



which has no 1/8 singularities. This form is similar to 
the value obtained in Ref. [28| using a rapidity regulator. 
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APPENDIX B: PARAMETER INTEGRALS 
1. Fermions 
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2. Scalars 
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For 4z > 1, the results can be analytically continued 
using Vl — 4z — ► iy/Az — 1 and tanh" 1 (Vl — 4z) — > 
i tan" 1 (V4z — 1) . In each integral, the factors of i cancel, 
and the function remains real. 
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